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Abstract 

We study dynamical and thermal effects that are induced in nanoparticle systems by a rotat- 
ing magnetic field. Using the deterministic Landau-Lifshitz equation and appropriate rotating 
coordinate systems, we derive the equations that characterize the steady-state precession of the 
nanoparticle magnetic moments and study a stability criterion for this type of motion. On this 
basis, we describe (i) the influence of the rotating field on the stability of the small-angle precession, 
(ii) the dynamical magnetization of nanoparticle systems, and (iii) the switching of the magnetic 
moments under the action of the rotating field. Using the backward Fokker-Planck equation, which 
corresponds to the stochastic Landau-Lifshitz equation, we develop a method for calculating the 
mean residence times that the driven magnetic moments dwell in the up and down states. Within 
this framework, the features of the induced magnetization and magnetic relaxation are elucidated. 

PACS numbers: 75.50.Tt, 76.20.+q, 05.40.-a 
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I. INTRODUCTION 



The study of the dynamics of the nanoparticle magnetic moments and their stability with 
respect to reorientations is a problem of prominent theoretical and practical importance. In 
fact, it is related to the stochastic and nonlinear dynamics and to the thermal stability of 
the magnetic moments in nanoparticle devices including magnetic storage ones?^ At low 
temperatures, when thermal fluctuations are negligible, the main interest is in the dynamics 
and stability in time-dependent external magnetic fields. In this case, the problem is usually 
reduced to the search for solutions of the deterministic Landau-Lifshitz equation 3 - and to the 
analysis of their stability. These investigations are also strongly motivated by the possibility 
of fast switching of the nanoparticle magnetic moments.- ' 5 ' 6 ' 7 ' 8 ' 9 

Due to thermal fluctuations, the dynamics of the nanoparticle magnetic moments be- 
comes stochastic and nonzero probabilities of their transition from one stable state to an- 
other appear. In this case, the dynamics can be described by the Fokker-Planck equation 
that corresponds to the stochastic Landau-Lifshitz equation.— At present this approach is 
widely used for studying magnetic properties of nanoparticle systems at finite temperatures, 
including magnetic relaxation , ^i 11 ' 12 ' 13 ' 14 ' 15 ' 16 

In this paper, we use the deterministic and the stochastic Landau-Lifshitz equations to 
study some effects induced by the rotating magnetic field in systems of purely deterministic 
and weakly superparamagnetic nanoparticles. More precisely, we are interested in the ef- 
fects that arise from the different dynamical states of the up and down magnetic moments. 
These states are generated by the magnetic field rotating in the plane perpendicular to the 
up-down axis and they are different even if the static magnetic field along this axis is absent. 
The reason is that the magnetic moments have a well-defined direction (counterclockwise) 
of the natural precession, and so the rotating field effectively interacts only with the up 
or down magnetic moments. We note in this context that some properties of the solutions 
of the deterministic Landau-Lifshitz equation were previously considered in the context of 
ferromagnetic resonance,— nonlinear magnetization dynamics , 18 ' 19 and switching of magne- 
tization in cylinders^ and spherical nanoparticles.— But, to the best of our knowledge, the 
above mentioned effects have not been investigated before. 

The paper is organized as follows. In Sec. II, we describe the model and the underlying 
assumptions. In Sec. Ill, we reduce the deterministic Landau-Lifshitz equation to the alge- 
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braic equations that describe the steady-state forced precession of the nanoparticle magnetic 
moments and derive a criterion of its stability. In the same section, we apply the results 
for studying the small-angle precession and switching of the magnetic moments. The effects 
in nanoparticle systems that arise from the simultaneous action of the thermal fluctuations 
and rotating field are considered in Sec. IV. Here we calculate the mean residence times 
that the driven magnetic moments reside in the up and down states, and apply these results 
to study the induced magnetization and magnetic relaxation in nanoparticle systems. We 
summarize our findings in Sec. V. Finally, in the Appendix we specify the used coordinate 
systems. 



II. DESCRIPTION OF THE MODEL 



We consider a uniaxial ferromagnetic nanoparticle with spatially uniform magnetization 
which is characterized by the anisotropy field H a and the magnetic moment m = m(i) 
of fixed length |m| = m. The assumption of uniform magnetization is valid for uniform 
nanoparticles if the exchange length, i.e., the length scale below which the exchange inter- 
action is predominant (for typical magnetic recording materials its order of magnitude is 
5-10 nm), exceeds the nanoparticle size. In other cases, e.g., for coated nanoparticles, it 
can be considered as a first approximation. We assume also that the static magnetic field 
H is applied along the easy axis of magnetization (the z axis), and the circularly polarized 
magnetic field h(t) is applied in the xy plane (see Fig. 1), i.e., H = He z and 

h(i) = h cos(u)t)e x + phsm(ujt)e y . (2-1) 

Here e x , e y , and e z are the unit vectors along the corresponding axes of the Cartesian 
coordinate system xyz, h = |h(t)|, uj is the frequency of rotation of h(i), and p = — 1 or 
+ 1 that corresponds to the clockwise or counterclockwise rotation of h(i), respectively. We 
write the magnetic energy of such a nanoparticle as 

W = -m 2 -Hm z -m-h(t), (2.2) 

where m z = m • e z is the z component of m, and the dot denotes the scalar product. 

In the deterministic case, we describe the dynamics of the nanoparticle magnetic moment 
by the Landau-Lifshitz equation 3 

rh = —7m x H e // m x (m x H e //). (2.3) 

777 
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Here 7(> 0) is the gyromagnetic ratio, A(> 0) is the dimensionless damping parameter, the 
cross denotes the vector product, and 

H e// = - — = h(t) + (H a ^ + H)e z (2.4) 
am V m J 

is the effective magnetic field acting on m. 

If the magnetic moment interacts with a heat bath, we use the stochastic Landau-Lifshitz 
equation^ 

rh = —7m x (Heff + n) m x (m x H e ff), (2.5) 

where n = n(t) is the thermal magnetic field with zero mean and correlation functions 
{n a {ti)np(t2)) = 2A5 a/ 35(t2 — £1), n a {t) (a = x,y,z) are the Cartesian components of n(t), 
A = Afc^T /7m is the intensity of the thermal field, is the Boltzmann constant, T is the 
absolute temperature, 5 a p is the Kronecker symbol, 5(t) is the Dirac 5 function, and the 
angular brackets denote averaging with respect to the sample paths of n(t). According to 
Eq. (|2.5jl . the conditional probability density P = P(6,ip,t\8' ,<p' ,t') (t > t') that describes 
the statistical properties of m in the terms of the polar 9 and azimuthal (p angles, satisfies 
the (forward) Fokker-Planck equatio n 1 ^ 22 

|p = _| (/l+f A cot() )P + 7 ^P 

-JL P (2.6) 
dip sin 9 o^p 2 

with 

m sin 9 \ o9 o<p I 



f2 = —^(ane^-X^-)w. 
m sin 9 \ 69 o<p J 



(2.7) 



III. DYNAMICAL EFFECTS 

A. Equations for the forced precession 

To study the forced precession of the nanoparticle magnetic moment and its stability with 
respect to small perturbations, we use the Landau-Lifshitz equation (|2.Hjl and represent m(t) 
in the form 

m(t)=m (t)+m 1 (t), (3.1) 



where m (t) describes the steady-state precession of m(t), and mi(t) is a small deviation 
from mo(t). Since |m(t)| = |mo(t)| = m, it is convenient to introduce the unit vector 
u = m (t)/m and a small dimensionless vector v = m 1 (t)/m (v = |v| <C 1). According to 
this, we decompose the effective magnetic field (|2.4j) into the zeroth-order (in v) vector 

= h(t) + {H a u z + H)e g (3.2) 

and the first-order one 

H« = H a v z e z . (3.3) 

Substituting Eq. (J3.1j) and the effective field H e ff = H-fff + h[jj into Eq. (|2.3|) and keeping 
the terms of the zeroth order, we end up with the following equation for u: 

li = - 7 u x H$, - A 7 u x (u x H$,). (3.4) 

Introducing, as usual, the rotating Cartesian coordinate system x'y'z' (see Fig. 1 and 
the Appendix) and assuming that in this coordinate system the components u x >, u y >, and 
u z i = u z of the vector u do not depend on time, Eq. (|3.4|) can be reduced to a system of 
algebraic equations. Indeed, using the relations 

u = u x >e x > + u y >e y , + u z e z , 

u = puj(-u y <e xl + u x >e y >), (3.5) 

= he* + (H a u z + H)e z 

that follow from Eqs. (jA.l|) - (|A.3|) and taking the x', y', and z components of Eq. ()3.4jl . we 
obtain 

\u x '{u 2 z + Hu z + hu x i) + u y i{u z + H — pu) = Xh, 

u x '(u z + H — pu) — \u y i{u 2 z + Hu z + hu x r) = hu z , (3-6) 

\u z {u 2 z + + hu x i) — \{u z + H) = huyi, 

where H = H/H a , h = h/H a , uj = u/u r , and uj t = jHa. A simple analysis of this system 

shows that u x > and u y > are readily expressed through u z . 

1 — u 2 Xk 
u x > = -= — -(uz + H-pn), Uyi = (1 u 2 ) (3.7) 

with k — uj/(1 + A 2 ), and u z satisfies the equation 

1 - v 2 

h 2 = -^±[( Uz + H - pn) 2 + (Xku z ) 2 ). (3.8) 

u% 



It is not difficult to verify that Eqs. (|3.7|) and ()3.8|) preserve the condition u 2 = 1. Note also 
that the components u x and u y of u in the initial coordinate system xyz are expressed in 
terms of u x > and u y i as follows: 

u x = u x > cos (cut) — puyi sin (tut), 

(3.9) 

u y = pu x i sin(u;t) + u y i cos(cjt). 

B. Stability criterion 

Next, assuming that the solution of Eq. (|3.8|) is known, we derive a stability criterion for 
the steady-state precession of m. To this end, using Eqs. (|2.3|) . (|3.4j) and (|3.1j) . we write the 
linear differential equation 

v = - 7 v x - 7 u x H ^ - A 7 [v(u • H^) 

+ u(v ■ H^) + u x (uxH«)] (3.10) 

that describes the evolution of small deviations v. Since Eq. (12. 3J) conserves |m|, the condi- 
tion 2u • v + v 2 = always holds. This means that, with linear accuracy in v, the vector v 
is perpendicular to u for all t. Therefore, it is convenient to introduce the rotating Carte- 
sian coordinate system x"y"z" (see Fig. 1 and Appendix) in which the vectors u and v are 
represented as 

u = e 2 » , v = v x » it) e x » + V y n it) e y n , (3.11) 

and so the condition u ■ v = holds automatically. 

In this coordinate system, the z" component of Eq. ()3.10|) is satisfied identically because 
according to the condition d(u ■ \)/dt = and Eq. (J3.4)) the relation 

[v + 7 v x H$/ + A 7 u(v • H^)] • e z „ = (3. 12) 

always takes place. Projecting Eq. (|3.1U|) onto the x" and y" axes and using Eqs. ()A.4|) - (jA.6|) . 
as well as the results of the previous section, we obtain after straightforward calculations 

i) x ii = —XuJiV x " — uj2V y ", 

(3.13) 

Vyll = 0JSV X " — \UJ4Vyl' } 
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where uo n = u r u n and 

Q 1 = u* + —[H-pK(l-u 2 z )], 

Uz 

u 2 = l + — [H - P k{\ + \ 2 u 2 z )} 



(3.14) 

Cu 3 = u 2 z + — [H - P k(1 + A 2 w 2 )], 

u z 

cD 4 = 1 + — [H-pK(l -u% 

Uz 

Thus, in the first-order approximation, the stability of the steady-state precession of 
the nanoparticle magnetic moment is defined by the stability of the stationary solution 
v x ii = Vyii = 0, or the fixed point (0, 0), of the system ()3.13|) . A complete solution of the last 
problem is well known (see, for example in Ref. [23]), and is based on the analysis of the 
roots 

A 



of the characteristic equation (5 + Aa>i)(5 + X004) + 012^3 = corresponding to this system. 
In particular, a criterion of the asymptotic stability of the forced precession has the form 
Re 84. < or 



X(Cji + 004) > Rey \ 2 (u)4 — uji) 2 — 4(^2^3 • (3.16) 

In the following we apply the above general results to study the precessional dynamics in 
the cases of small precession angles and zero static magnetic field. 



C. Small precession angles 

In this case, we assume that the precession angles 9 a (see Fig. 2) of the magnetic moments 
with u z > (er = +1) and u z < (a = —1) are small, i.e., 6 a 1. Then u z can be 
represented in the form u z = er(l — e 2 /2), where according to Eq. (|3.8J) a small parameter e 2 
is given by 

h 2 

e 2 = ^ , (3.17) 

(a + H - pk) 2 + A 2 k 2 

and so Eq. ()3.7j) yields with linear accuracy in e 

1 + <jH — apK 



u r j = h- 



(a + H - P k) 2 + \ 2 k 21 

y H ! (3.18) 

~ p\n 

u„i = —a = . 

v (a + H - P k) 2 + \ 2 k 2 



We emphasize that, even though the static magnetic field is absent, i.e., H = 0, the dynamics 
of the up (cr = +1) and down (cr = —1) magnetic moments is quite different. The reason 
is that the natural precession of the magnetic moments is counterclockwise, and so only up 
or down magnetic moments have the direction of the natural precession that coincides with 
the direction of the magnetic field rotation. In other words, the magnetic field rotating in 
the plane perpendicular to the easy axis of magnetization breaks the degeneracy between 
the up and down states of the magnetic moment. 

Our analysis shows that the small-angle precession is stable only if £j\ + tl> 4 > 0. Writing 
Cj\ + with quadratic accuracy in e, 

Q 1 + u) 4 = 2 + 2aH - e 2 (l - oH + 2apn), (3.19) 

and solving the equation Q\ + cj 4 = with respect to H, we find the critical magnetic field 

that separates the stable and unstable precession for a given state a. The steady precession 
is stable either for H > H cr \ a=+1 or for H < if cr | (T= _ 1 . Because at H = H cr the precession 
in the state —a is stable, the switching of the nanoparticle magnetic moments from the 
unstable state a to the stable state —a occurs. If p = cr, then the rotating field always 
decreases the stability of the precession, i.e., H cr \ a=+1 > —1 and H cr \ a =-i < +1. On the 
contrary, if p = —a then depending on the reduced frequency k the rotating field can both 
decrease (if k < 1) and increase (if k > 1) the stability. The largest stabilization effect is 
achieved at k — 2. Note also that since e«l and usually A < 1, the formula (|3.20|) is valid 
only if the condition k ^> h holds. 

As an important illustrative example, we consider the nanoparticle system with the same 
number N/2 of the up and down states. In this case the dynamical (dimensionless) mag- 
netization of the system pd = (l/N)J2iLi u zi {i labels the nanoparticles) takes the form 
Pd = (1/2) J2a ° cos ®(T or ) since 9 a <C 1, p.d — (# 2 i _ 0+i)/4. Assuming that H = and 
using the formula 

e a = = (3.2i) 

^1 — apn) 2 + X 2 k 2 



that follows from the relation sin# CT = e and Eq. (j3.17j) . this quantity can be written in the 
form 

^ d = ~S + (l + A>2p-4^ (3 - 22) 
8 



This result shows that (i) the magnetization pd is a purely dynamical effect, i.e., pd = 
if k — 0/(1 + A 2 ) = 0, (ii) the direction of magnetization and the direction of magnetic 
field rotation follow the left-hand rule, and (iii) the dependence of pd on k always exhibits 
a resonant character. The maximum of fid occurs at k = K m , where 

Km = ^ + \A - A 2 + 2v / T+A 2 "T^, (3.23) 

and fi d \ K=Km = —p(h/2X) 2 for A < 1. If h -C A then the dynamical magnetization is small 
but, as we will show later, it can be considerably enhanced by thermal fluctuations. 



D. Zero static magnetic field 

In the case of zero static magnetic field, H = 0, we rewrite Eq. ()3.8|) in the form h = 
F p (u z ), where 

F„(x) = V ^ X ^/(x- P k) 2 + (\kx) 2 (3.24) 

(— 1 < x < 1). According to the definition, the function F p (x) satisfies the conditions 
Fp(—x) = F- p (x), F p (—1) = F p (+1) = 0, F p (x) — > oo as x — > 0, and it has a local minimum 
at x — px\ {x\ > 0) and a local maximum at x — px 2 (x 2 > 0), see Fig. 3. 

A detailed analysis shows that for fixed p the precession of the nanoparticle magnetic 
moment in the state a = —p is stable for all values of h. In other words, the unique solution 
of the equation h = F p (u z ) with sgnu 2 = — p always exists and is stable. In this case, the 
dependence of u z on h is shown in Fig. 4, curve 1. 

The precession of the magnetic moment in the state a = p (when sgnn 2 = p) exhibits 
qualitatively different behavior depending on h. It is stable only if \u z \ > x 2 that implies 
that h < h cr = Fp(px 2 ) (see Fig. 4, curve 2). At h = h cr + the solutions u z \ a=p = px 2 
and u z \ a=p = px% (x% > 0) are unstable, and the magnetic moment switches from the state 
with u z \ u=p = px 2 to the new state with u z \ a= ^ p = — px 4 (x 4 > 0). As stated above, the new 
state is stable for all h, and so the reverse transition does not occur for fixed p. 

In the important case of small driven frequency, k <C 1, that is easily accessible to 
experimental investigation, the analysis of the stability of the forced precession can be done 
analytically. In particular, we found that 



if 1 - h 2 > k 2/3 and u 



p ~ 



pn/h if h ^> 1. Using the approximate representations 



UJl = UJ 3 = III - p , U 2 = Co> 4 = 1 - p 

u z u z 



(3.26) 



we showed explicitly that the stability criterion ({SHE} for a = p is reduced to \u z \ > x 2 and 
the critical amplitude of the rotating field is given by h cr = 1 - (3/2)k 2 / 3 . Finally, solving 
the equations dF p (x)/dx = and h cr = F p (x) with respect to x, we derived the asymptotic 
formulas X\ = k, x 2 = k 1//3 , X3 = k/2, and X4 = 2k 1 / 3 . It is important to note that although 
the jump Au z = x 2 + £4 that occurs under switching of u z tends to zero as n — >• 0, it can be 
sizeable even for very small k (for example, Au z = 0.3 if k = 10 _3 ). 

We emphasize also that this truly remarkable phenomenon, the switching of the nanopar- 
ticle magnetic moments under the action of the rotating field, results from the existence of 
the natural precession of the magnetic moments. It occurs only in those nanoparticles for 
which the condition op = +1 holds, i.e., if the direction of the magnetic field rotation 
coincides with the direction of the natural precession of the magnetic moments. 

IV. THERMAL EFFECTS 
A. Mean residence times 

If the magnetic moments interact with a heat bath then their dynamics becomes stochastic 
and is described by the forward Fokker-Planck equation ()2.6|) . In this case, due to thermal 
fluctuations, the magnetic moments can perform random transitions from the one state a 
to the other —a. Our aim is to study how the rotating magnetic field influences the mean 
residence times t a that the magnetic moments dwell in these states at H = 0. In principle, 
the problem can be solved on the basis of Eq. (|2.fij) . Specifically, this approach has been used 
in the case of ac magnetic field linearly polarized along the easy axis of magnetization.— 
However, since the mean residence times can be readily expressed through the mean-first 
passage times, for solving this problem it is convenient to use the backward Fokker-Planck 
equation22^ 




(4.1) 
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[fi 2 = f 1,2(6', (f',t')}, which is equivalent to Eq. (|2.6|) . Within this framework, we are able 
to calculate t a in some particular cases. But because of the procedure is rather complicated 
from the mathematical point of view (details will be published elsewhere), here we use a 
crude approximation that leads, however, to qualitatively the same results. 

In the considered case of small-angle precession, the magnetic moments of weakly super- 
paramagnetic particles (when a = H a m/2kBT 3> 1) spend almost all time near the conic 
surfaces with the cone angles (|3.2ip . We assume that if these imaginary surfaces are replaced 
by the reflecting surfaces, then the rotating field terms can be eliminated from Eq. (|4.1jl . 
Then, replacing also the conditional probability density P by P = P(8,t\6' ,t') and taking 
into account that f[ = — (A/2)u; r sin 2$', Eq. (|4.1|) reduces to the simpler form 

P = (* . r sin 20' - 7 2 A cot f a|>. (4.2, 

Next we use a standard procedure 2 ^ to define the mean first passage times for the magnetic 
moments in the states a 

roc r e^ 

T a (6') = / dt d6P(6,t\6',0) (4.3) 

and to derive from Eq. f|4.2jl the ordinary differential equation for these quantities 

+ (cot? - a sin 29') ^1 + at r = 0. (4.4) 
a9' A a9' 

Here, & G (9 +1 , tt/2) if a = +1, & G (vr/2, tt - 9_ x ) if a = -1, 0$ = O+i, 9 +l = e -i = 

9_\ = tv — 6-i, and t r = 2/\u r is the characteristic relaxation time of the precessional motion 

of the magnetic moment. Solving Eq. (|4.4|) with the absorbing and reflecting boundary 

conditions, i.e., T a (n/2) = and dT a (6') / d6'\ e i =7l (i_ a y2+ae a = 0, respectively, and taking 

into account that the desired times t a are readily expressed through T a (6'), t a = 2T a (n[l — 

cr]/2 + a6 a ), we obtain for a ^> 1 and 6 a ^ 1 

t a = exp[a(l + o2H ff )\, (4.5) 

where H a = —o6 2 a j2 can by interpreted as an effective magnetic field acting on the nanopar- 
ticles in the state a. 

According to this result, the rotating magnetic field decreases the mean residence times. 
However, due to the natural precession, the decrease of the mean times is different for the 
up and down magnetic moments. As it will be shown below, this fact causes a strong 
enhancement of the dynamical magnetization and leads to a modification of the relaxation 
law. 
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B. Induced magnetization 



We define the steady state magnetization of the nanoparticle system in the rotating 
magnetic field as /i — (1/iV) Si=i Tn z ijm (N — > oo). Denoting the average number of 
the magnetic moments in the state a as N a and introducing the probability p a = N a /N 
(p+i + P-i — 1) that the magnetic moment resides in this state, we rewrite /i in the form 

/i = 5> CT (cos0>| CT , (4.6) 

a 

where (cos^)^ = (1/N a ) J2i^N a cos ®i is the average value of m zi /m in the state a. If a 3> 1 
then thermal fluctuations are small and (cos#)| CT in Eq. (J4.6)) can be replaced by a cos Q„ ~ 
cr(l - 9l/2), yielding fi = fH + fkd + V>d- Here fi t = E CT °Vo and m & = £ CT o"(2p CT - 1)6^/4 are 
the contributions of thermal fluctuations to the total magnetization /i, and /i<j is the purely 
dynamical magnetization given by Eq. (|3.22j) . Since 6 a <C 1, the condition m] [ltd ^ 1 
holds, and so \i ~ \i t + /i^. We emphasize that a decrease of the temperature (increasing 
of a) decreases the fluctuations of the magnetic moments, but not the difference between 
p + i and i.e., // t . Moreover, one expects that, similar to the two-level models, |// t ] grows 
with a. 

Next, taking into account that in the steady state p a = t a j and using Eq. ()4.5|) . 
we obtain 

Ht = tanh[a(# +1 + H-i)}. (4.7) 

Comparing fi t with the magnetization of an Ising paramagnet, tanh(2aif), we see that the 
circularly polarized magnetic field induces the same magnetization \i t of the nanoparticle 
system as the external magnetic field (if + i+if_i)/2 applied perpendicular to the polarization 
plane. It is interesting to note that (H+i + H-i)/2 = fid, and hence /i t = tanh(2a/Xd)- Since 
a ^> 1 and \fid\ <C 1, this relation shows that fa/l^d ^ 1 an d so u ~ /if, i.e., thermal 
fluctuations strongly enhance the dynamical magnetization. In particular, if a\fid\ ^ 1 then 
Ht/l^d = 2a. Like the dependence of /if on « has a resonant character and, as expected, 
| /it | increases with decreasing temperature (see Fig. 5). 

C. Relaxation law 

As a second example, let us consider the thermally activated magnetic relaxation in the 
nanoparticle system driven by the rotating field. Since the transition rate of the nanoparticle 
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magnetic moment from the state a to the state —a equals l/t a , the differential equation that 
defines the time-dependent magnetization fi{t) of this system can be written in the form 

fl(t) = -fi(t) ( 7^ + rH " 7^ + 7^- (4-8) 

\t+i t_i j t+i r_i 
Assuming that /i(0) = 1 (we neglect the dynamical magnetization), from Eq. ()4.8|) we obtain 
the relaxation law 

fi(t) = (1 - fi t ) exp(-t/r) + fi u (4.9) 

where 

r= expK^-if,,)] 
cosh[a(F +1 + iZLj)] 

is the relaxation time in the presence of the rotating magnetic field, and r = 



(t r /4)J ir/aexp a is the relaxation time if the rotating field is absent. 

Thus, the rotating magnetic field decreases the relaxation time (t/t < 1) and leads to 
nonzero magnetization in the long-time limit [fi(oo) = fit]. 



V. CONCLUSIONS 



We have investigated a number of dynamical and thermal effects in nanoparticle systems 
that result from the action of a circularly polarized magnetic field rotating in the plane 
perpendicular to the easy axes of the nanoparticles. The main finding is that the dynamics 
of the nanoparticle magnetic moments, both deterministic and stochastic, becomes different 
in the up and down states. It is important to note that, due to the (counterclockwise) 
natural precession of the magnetic moments, the dynamics is different even if the static 
magnetic field is absent. 

To describe the dynamical effects at zero temperature, we have used the deterministic 
Landau-Lifshitz equation. We have solved this equation for small-angle precession of the 
magnetic moments and have demonstrated that the rotating field, depending on its frequency 
and polarization, can either decrease or increase the stability of the precession motion. For 
zero static field, we have calculated the dynamical magnetization of nanoparticle systems 
and predicted the switching effect. This remarkable effect, which consists in changing the 
state of the magnetic moments at some critical amplitude of the rotating magnetic field, 
occurs only for resonant nanoparticles, i.e., when the direction of the natural precession of 
their magnetic moments coincides with the direction of the magnetic field rotation. 
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In the case of finite temperatures, we have invoked the backward Fokker-Planck equation 
to calculate the mean residence times that the driven magnetic moments dwell in the up and 
down states, respectively. On this basis, we have studied the steady-state magnetization and 
the features of magnetic relaxation in systems of weakly superparamagnetic nanoparticles 
that are driven by the rotating magnetic field. In particular, we have found that thermal 
fluctuations strongly enhance the dynamical magnetization and that the rotating field always 
causes a decrease of the relaxation time. 
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APPENDIX: ROTATING COORDINATE SYSTEMS 
1. Single-primed coordinate system 

The rotating Cartesian coordinate system x'y'z' is defined by the unit vectors e x >(t), e y i(t), 
and e z > (t) that are expressed through the unit vectors of the initial (laboratory) coordinate 
system xyz as follows: 

e x i = cos(ut)e x + p sin {ujt) e y , 

(A.l) 

B y i = —psin(u>t)e x + cos(ut)e y , 
£z' = e 2 . According to Eqs. (jA.l|) . the inverse transformation has the form 

= cos(u;t)e x / — psm(ut)e y i, 

(A.2) 

e y = p sm(idt)e x > + cos (ut) e y >, 

and 

e x > = pue y/ , ey = -pue x >. (A. 3) 
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2. Double-primed coordinate system 



The unit vectors e x ii(t), e y n(t), and e z »(t) of the rotating Cartesian coordinate system 
x"y"z" are introduced as 



Vx" = , „ = Ju 2 x , +u 2 ,e z ,, 



V U x' + u y' 



u y ie x i + u x ie y i ^ ^ 



e z " = u x >e x i + UyiQyi + u z 'B z i, 
and so the inverse transformation is given by 

U X 'U Z /6 X " UylGyll 

&X 1 — ; 1" Ux'^z", 

<< + s> 



u y >u z >e x r> + u x >e y >i i . -x 

<V = —T= — + Uy>e z „ } (A.5) 



^u 2 x , + u 2 , 

e z i = -\ju\, + tty/ ea,« + w z 'e 2 //. 
From here and (jA.3|) . straightforward calculations yield 

e x ii = pUJUz'Gy", 



e y ii = —poo iu z >e x » + Ju 2 , + u 2 , e z n 



(A.6) 
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FIG. 1: Schematic representation of the model and the used coordinate systems. 
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FIG. 2: Sketch of the precession angles for the up and down magnetic moments (the arrows depict 
the directions of their natural precession). 
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FIG. 3: Plot of the function F +1 (x) for k = 0.25 and A = 0.2. If k -> 0, then x n -> 0, F +1 {x 1 ) -> 0, 
and h cr = F+ifaz) — > 1. 
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FIG. 4: Plots of the stable solutions of the equation h = F+i(u z ) for the same parameters as in 
Fig. 3. For p = — 1 the curves 1 and 2 must be reflected with respect to the axis h. 
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FIG. 5: Plots of the magnetization fit vs the reduced frequency k for the parameter choice: 
h = 10 -3 , A = 10~ 2 , a = 20 (curve 1) and for a = 40 (curve 2). The temperature in the latter case 
is two times less than in the former one. 
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